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Abstract. We develop a classical theoretical description for nonlinear many-body 
dynamics that incorporates the back-action of a continuous measurement process. The 
classical approach is compared with the exact quantum solution in an example with an 
atomic Bose-Einstein condensate in a double-well potential where the atom numbers 
in both potential wells are monitored by light scattering. In the classical description 
the back-action of the measurements appears as diffusion of the relative phase of the 
condensates on each side of the trap. When the measurements are frequent enough to 
resolve the system dynamics, the system behaves classically. This happens even deep 
in the quantum regime, and demonstrates how classical physics emerges from quantum 
mechanics as a result of measurement back-action. 
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It is now possible to carry out experiments with systems comprised of a small 
and controllable number of atoms, and eventually the interface between quantum 
mechanics and classical physics [H El E] must be confronted. Measurement back- 
action is an essential element of quantum mechanics, and may play a key role also 
in the transition from quantum mechanics to classical physics |11 El [6], |7]. In this 
paper we introduce a theoretical method to account for measurement back-action in 
nonlinear dynamics [H |9l [TOl [HI [I2l US] classically, using a continuously monitored [H] 
Bose-Einstein condensate (BEC) in a double-well trap as an explicit example. When 
the continuous measurements are strong enough so that the system dynamics can be 
resolved, the classical description and a full quantum analysis agree well even deep 
within what might appear to be the domain of quantum mechanics. 

We employ the minimal two-mode model [15] to describe the Josephson oscillations 
of the atoms between the potential wells [161 113 HI [El 120] , see figure [H Given the 
hopping amplitude J between the wells, the strength of atom-atom interaction inside 
a well U, and the imbalance in the energies of the atoms between the sides e, the 
Hamiltonian is 

H = -J{a^b + b^a) + U{a^a^aa + b'^b^bb) + e{a^a - b^b) . (1) 

Here a and b are the boson operators for the atoms in the two wells, and we set /i = 1. 

The atom numbers in the two wells are monitored by coherent scattering of off- 
resonant light; consult figure [H Counting of photons on two detectors, one for each 
potential well, constitutes the final measurement. In our quantum description we 
consider the atoms as an open system, and model the back-action of the measurements 




Figure 1. Left: Bose-Einstein condensate in a double-well potential with tunneling 
amplitude J, atom-atom interaction strength C/, and energy difference per atom 
between the wells e indicated. Right: Basic idea of detection of the atoms in one 
trap. Since the detection light is far-off resonance, the dominant effect of the atoms is 
to alter the phase of the light passing through the sample by an amount proportional 
to atom number. The phase shift is then detected interferometrically. The amplitude 
of the light passing to the photodetector is proportional to the number of atoms, so 
that the rate of photon counts is proportional to the square of the atom number. The 
scheme is heavily simplified. For instance, mode matching of the light to the condensate 
is not considered. 
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by adding a linear relaxation term to the equation of motion of the density operator of 
the atoms, 

p = -z[H,p] + Cp. (2) 

The most general relaxation term Cp such that a density operator remains valid during 
the time evolution is of the Lindblad form [21], 122] . 

Cp = J2 [l^p4 - UpLlL^ + LlL.p)] , (3) 

i 

where Li are operators of the double-well system. In the present model [H] we have the 
Lindblad operators La = y/V a^a and Lb = y/Tb^b, indicating that the scattering rates 
of photons from the two wells are T{{a^a)^) and T{{b^b)^), respectively. If the light is 
predominantly scattered into a narrow cone in the forward direction, F also specifies the 
strength of the measurements. It can be adjusted, for instance, by varying the intensity 
or the frequency of the probing light. 

One could develop a parallel classical description using various phase-space 
representations such as the Wigner function for the ordinary [231 1211 |22] or the 
SU(2) [251 EEl ET] coherent states, or the number-phase Wigner representation [28] . 
We choose the standard Wigner function because it is generally used in quantum optics 
and is not specifically tailored to the present system. We begin by rewriting the master 
equation ([2]) equivalently in terms of the Wigner function W{a,a*, (3, f3*;t) [211 122] , 
where a and (3 are complex numbers corresponding to the operators a and b. In the 
equation of motion of the Wigner function, a, a* and (3, /3* behave like canonical pairs 
with Poisson brackets such as {a, a*} = —i. We then change to four real variables N, 
$, Z and if defined by 



a = ^ \N - Z e'^'^-~2^\ 13 = ^ \N + Z e''-'^+2^^ . (4) 

This transformation between the variables is canonical, and results in Poisson brackets 
such as {Z, if} = 1. To the extent possible, we regard the Wigner function as the classical 
joint probability density for total atom number A^, imbalance of atom numbers between 
the sides of the double well Z {Z E [—N/2, N/2]), overall phase of the two BECs $, and 
phase difference of the BECs between the two potential wells ip. The global phase $ 
may be dropped as it has no dynamics, nor any relevance to the dynamics of the other 
variables. We define x = UN, z = 2Z/N {z G [—1,1]), and expand the equation of 
motion of the Wigner function up to and including the order 0{N^^) in the formally 
large atom number N. The equation of motion 

9h g ^ dh 
dip ^ dz 

follows for the remaining distribution function W{z,(p), with the effective one-particle 
Hamiltonian 



dtW 



W (5) 



h{z, (p) = z^x - 2ez - 2 JVl - z"^ cos((^) . (6) 

Importantly, the quantum mechanical back-action of the measurements of the 
particle numbers is now accounted for by classical diffusion of the relative phase of the 
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BECs proportional to Td^^. This is, perhaps, as expected: The relevant measurement 
result is the difference of the photon counting rates on the two detectors, and this is 
proportional to the expectation value of the operator 

AN = {a^af - {h^hf = {a^a + b^b){a^a - b^b) . (7) 

But the first factor in the final form is the conserved total atom number, so the measured 
operator is actually the difference of atom numbers between the sides of the trap. Given 
the (in quantum mechanics, qualitative) conjugate relation between number and phase, 
measurement of atom number difference should indeed increase the uncertainty of the 
relative phase. 

We compare the quantum and classical evolutions in the case when the system is 
first prepared to the ground state at zero temperature in the presence of an energy 
imbalance, e 7^ 0, and the imbalance is suddenly removed by setting e = 0. The atoms 
then start oscillating between the two sides of the potential well. 

In our quantum analysis the evolution of the density operator according to ([2]) is 
solved using quantum trajectory simulations [29l [301 EI] . For every run of the simulation 
we have a state vector that undergoes 'quantum jumps' at random times according 
to a precisely defined probabilistic law. As discussed in, say, Refs. |22l [32], here we 
also assume an "efficient" measurement process in which every scattered photon gets 
detected. Then a quantum trajectory simulation with its quantum jumps produces a 
faithful simulation of an individual experimental run, each quantum jump corresponding 
to a detection of a photon. Even though the master equation ([2]) is not conditioned on 
the results of the measurements and as such describes an experiment in which the 
measurement results are discarded, the results are built into each quantum trajectory. 

In a real experiment one would record a stream of photon counts, and the population 
imbalance z{t) has to be inferred from them in some manner. For instance, if frequent 
photon counts merge into an essentially continuous photocurrent that varies on a time 
scale of the back-and-forth oscillations of the atoms, it would be sensible to regard 
the difference of the photocurrents coming from the two detectors as a measure of 
the instantaneous population imbalance. However, for the sake of a clean illustration 
we analyze our quantum simulations differently: Given a quantum trajectory, a time 
dependent state vector, we use the expectations value of z = {a'^a — b'^b)/N to represent 
the value of the experimentally inferred population imbalance z{t). Evidently this works 
the better, the more numerous the quantum jumps are on the time scale of the variation 
oizit). 

In the classical analysis we take into account the quantum fluctuations in the 
initial state in the same way as is customary with the truncated Wigner approximation 
(TWA) [331 1311 [35l |36] . Thus, suppose equations (jS]) and (j6]) were written in terms of the 
variables Z and instead of z and if. We expand the corresponding Hamiltonian ^{Z, ip) 
to second order in Z and if around the classical lowest-energy state in the presence of the 
initial energy imbalance e 7^ between the sides. Since the classical variables Z and 
are a canonical pair, we may requantize by postulating the commutator [Z, ip>] = i. 
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If) thereby becomes the quantum Hamiltonian of a harmonic oscillator. The 
Wigner function corresponding to the ground state of this oscillator is used as the 
initial distribution of the variables Z and ip, and leads to the initital condition for the 
distribution of the variables z and (f 



N 

Wo{z,ip) = — exp 



(8) 



2 V i 

We take x > and choose e in such a way that the classical equilibrium has the phase 
= and a predetermined population imbalance z = f. The parameter 

1 - P 

^ = (9) 

+ (1 - pr^'x/j 

reflects the size of the ground state wave function of the underlying harmonic oscillator. 

The evolution with the initial Wigner function Wo{z, ip) is solved by unravelling 
the dynamics into stochastic trajectories [33l [311 [371 [26] . If there were no phase diffusion 
in ([2D, we would have an ensemble of trajectories {z{t), v?(t)) evolving under the effective 
Hamiltonian h, which makes the usual TWA. Given the diffusion, there exists a precise 
mathematical correspondence between and the stochastic differential equations 



dz = 2JV1 - z"^ sin dt, (10) 
dip= -2z(x + J . ^"^'^ ) dt + V2fdW{t) , (11) 



where dW{t) denotes a random increment of the Wiener process over the interval dt |37]. 
The evolution of the distribution W{z, ip; t) from the initial condition VFo(^? v) according 
to (12]) is obtained from an ensemble of stochastic trajectories. Specifically, the initial 
state of each trajectory {z,ip) is sampled from the distribution Wq{z,ip), the dynamics 
of the trajectory {z{t), (p{t)) is determined from f[TO]) and f [TT]) . and the expectation value 
of any quantity over the distribution W{z, ip; t) is correctly obtained by averaging the 
value of the quantity over a large number of such trajectories. The stochastic dynamics 
is integrated numerically using the Milstein algorithm |37] . 

We are now in a position to state precisely what "classical description" means in 
this context: unlike quantum mechanics, the analysis conforms to classical logic. The 
usual mean-field theory [T7] is a classical description, as are the TWA and our present 
approach. The initial Wigner function Wq is a valid classical probability distribution, 
even when it approximately incorporates quantum fluctuations of the atomic populations 
in the initial state. Moreover, subsequent evolution according to ([5]) that also accounts 
for measurement back-action keeps the probability distribution function valid at all 
future times. Correspondingly, we have been able to represent the complete solution 
in terms of classical stochastic trajectories. Such a trajectory exists objectively even 
without observations, but otherwise a similar correspondence between simulations 
and experiments prevails in the classical description as in quantum mechanics: each 
individual stochastic realization would be a representative outcome of an individual 
experiment. 
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Figure 2. An individual classical (upper trace, blue) and quantum (lower trace, red) 
trajectory. The parameters are N = 100, x/^ = 10, / = 0.5, and F = 0.003. 

Figure |2] shows a classical (blue upper trace) and a quantum rendition (red lower 
trace) of a single experiment. In the mean-field theory, and for each individual TWA 
trajectory, one expects strictly periodic oscillations in the classical trace. The deviations 
from periodicity originate from the phase diffusion. On the other hand, the measurement 
strength V / J = 0.003 is chosen such that each photon detector produces about ten clicks 
during a typical back-and-forth swing of the atoms, so that the quantum trace is a fair 
representation of a measurement record. We have deliberately picked these random 
examples so that qualitatively similar classical and quantum behavior with occasional 
trapping of the populations |17] in one or the other well is seen. 

Nevertheless, even if we are able to model individual experiments both classically 
and quantum mechanically, these are stochastic stimulations. In order to construct 
expectation values of experimental observables for quantitative comparisons we have 
to average over many trajectories. In figure [3] we demonstrate simulation results for 
population imbalance when averaged over a thousand trajectories. The upper two panels 
are for = 100 atoms and the lower two for = 10, while the differences between the 
graphs (a)-(d) are in the measurement strength F. 

Panel (a) for A^ = 100 atoms plots the quantum and classical expectation values 
of the population imbalance z{t) with no measurements, F = 0. The quantum result 
then represents the usual quantum average of the operator z = {a'^a — b'^b)/N, and 
has nothing to do with the (nonexistent) continuous monitoring. The quantum and 
classical graphs both show a collapse up to Jt ~ 10, but only the quantum expectation 
value has a revival at around Jt = 70. In panel (b) we turn up the measurements 
strength to F/J = 3 x 10^^. Both detectors then report about ten photons during a 
characteristic oscillation time of the atoms between the traps, so that the oscillations 
of the atoms could be resolved by monitoring the light scattered off the atoms. The 
quantum expectation value of z should therefore be a passable representation of the 
population imbalance inferred from photon counts. The collapse part, quantum and 
classical, is almost unchanged by the measurements, but the quantum revival has 
disappeared. The quantum and classical analyses now agree very well. 

We finally move on to the atoms number A^ = 10. In panel (c) we show the 
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Figure 3. Quantum (solid red lines) and classical (dashed blue lines) expectation 
values of the population imbalance 2: as a function of time t, all averaged over 1,000 
trajectories. The upper half of the figure is for atom number N = 100, the lower half 
for N = 10. Graphs (a) and (c) are without continuous monitoring, F/J — 0, graphs 
(b) and (d) correspond to the observation strengths T / J — 0.003 and F/J = 0.3 . The 
atom-atom interaction strength is xl J = lOi ^-nd the initial population imbalance is 
/ = 0.5. The insets demonstrate the quality of the quantum-classical agreement in the 
collapse region of traces (b) and (d) . 

population imbalances as a function of time for both the quantum theory and the 
classical theory without continuous monitoring, F = 0. There is little in common 
with the two traces, and one would conclude that the classical description is totally 
off. The proximate reason is that with a decreasing atom number the quantum revivals 
start overlapping and produce irregular-looking behavior that is absent in the classical 
approach. It appears that we are deep in the domain of quantum mechanics. However, 
in panel (d) we have continuous measurements with the strength T/J = 0.3, which 
again gives about ten photon counts during a back-and forth oscillation of the atoms; 
and again, the classical and quantum results agree. 

We have surveyed the behavior of the system for parameter values comparable 
to those in figures [2] and El In the absence of measurements, F = 0, we have 
the usual TWA and the corresponding full quantum theory. The TWA relies on 
a stochastic description and cannot fully incorporate quantum effects. The generic 
observation is that TWA agrees well with the exact solution for short times when 
stochastic field is highly occupied, but becomes progressively less accurate a model 
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of quantum dynamics when the average mode population becomes smaller and the 
simulation time scales longer [331 EH |35] . There are also examples, however, especially 
in lower dimensional systems in which cases the TWA simulations have been successful 
in providing qualitative descriptions of experimental findings even in the strongly 
fluctuating quantum limit [38] . 

In figure |3], in traces (a) with = 100 and F = we find a reasonable agreement 
between the TWA and quantum solutions up to the quantum revival around Jt = 70, 
while in figure |3]^c) with A^ = 10 the two cases part ways after one back-and forth swing 
of the atoms at Jt ~ 1, a time that is hard to discern in the figure at all. Second, as 
the measurement strength is increased from zero, the changes in the quantum results 
start being qualitatively visible at times when there have been, on the average, a few 
counts on both detectors. The classical results also change with increasing F, but the 
changes are less conspicuous. For example, with A^ = 100 the maximum amplitude 
of the quantum revival at Jt ~ 70 has decreased to half of its value in figure [3]^a) for 
F/J = 2 X 10~^, at which point the average time between photon counts at each detector 
is about 20/ J. Let us next consider F/J = 3 x 10^'^, the value used in figure [3]^b) for 
A^ = 100. Since our definition of the parameter F is such that the photon counting rate 
on the detectors is proportional to the square of the atom number, for A^ = 10 this 
measurement strength is still modest, with the time between counts of about 10/ J. A 
plot (not shown) would reveal that by Jt = 100 the oscillations in the quantum result 
have lost 2/3 of their amplitude, but that the quantum and classical approaches still 
give quite different results. 

Finally, as soon as the atom numbers are monitored vigorously enough so that the 
classical evolution in the form of the oscillations of the atoms between the two sides of the 
double- well trap can be resolved in detail in the scattered light, figures [3]^b) and (d), the 
quantum and classical descriptions give materially the same results. In round numbers, 
the time scale for classical evolution in our examples is 1/ J and the scattering rate to 
both detectors is A^^F/4. The watershed criterion is therefore N'^T/AJ ~ 1, exceeded 
by a factor of about 10 in both graphs [3](b) and (d). Moreover, and unlike in the 
usual case with the ordinary TWA, in our amended version of the TWA the quantum- 
classical agreement seems to persist for a long time, possibly as long as our numerical 
computations remain reliable. 

The emergent classicality is a consequence of measurement back-action suppressing 
the subtle interference effects associated with quantum revivals that are inherent in 
the full quantum description [3]. Remarkably, this occurs even deep in the domain of 
quantum mechanics with only 10 atoms, representing a far more dramatic transition to 
classicality than, e.g., measurement-induced relative phase between condensates [5l[6l[7] 
or nonlinear instabilities [TOl [121 1131 EH]- The classical description at this level of 
accuracy is only possible because we have included classical phase diffusion to account 
for quantum mechanical measurement back-action. To wit, the only difference between 
the classical predictions in graphs M^c) and (d) is the amount of phase diffusion, yet 
they are clearly different. 
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The main technical achievement of this paper is that we have introduced a 
theoretical method for a continuously monitored quantum system that handles both 
quantum fluctuations and measurement back-action classically. Now, the computational 
resources required to solve the time evolution from quantum mechanics tend to increase 
very fast with the system size. We therefore hope to inspire practical algorithms to 
analyze time evolution in situations where quantum effects are small but not negligible. 
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